Vortex-induced strain and magnetization in type-II superconductors. 
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It is argued that the stress caused by vortex cores in the mixed state of superconductors may 
result in a field dependent contribution to the free energy and magnetization of measurable levels. 
For sufHciently strong stress dependence of the critical temperature, dT^/dp, this contribution may 
result in the so-called "second peak" in the field dependence of the reversible magnetization, the 
effect often masked by vortex pinning and creep. 
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The so-called second peak in the field dependence of 
the magnetization M{H) in a number of type-II super- 
conductors is a long standing puzzle. The peak has 
been observed in magnetization loops having a "fish tail" 
shape so that the loop width of the irreversible magneti- 
zation increases with increasing field in intermediate field 
range suggesting the critical current rise with increasing 
field in this domain.^ 

A few explanations based on peculiarities of pinning 
and flux creep have been offered for this apparently 
strange phenomenon.!^ These suggestions may well be 
correct but they do not cover all cases in which the sec- 
ond peak has been observed. Puzzling in particular is 
the fact that in some systems the second peak has been 
reported also in the reversible M{H); the exaniples are 
NbSes, Lai.45Ndo.4oSro.i5Cu04, and CeCoIng.l^ 

In this work, the second peak in reversible M{H) is 
associated with the strain caused by normal vortex cores 
embedded in the superconducting phase, a "magneto- 
elastic" effect. The strains arise due to a small difference 
in densities of the normal and superconducting phases 
which is related to the stress dependence of the crit- 
ical temperature dTc/dpW It turned out recently that 
this derivative in pnictides, and in Ca(Fei_3;COa;)2As2 in 
particular,'^ by one or two orders of magnitude exceeds 
values for conventional superconductors making Fe-based 
pnictides especially favorable for observation of magneto- 
elastic effects. 

Strain caused by a single vortex. Consider vortex 
nucleation prior to which the superconductor has been 
strain free. We model the vortex core as a normal (n) 
cylinder of radius p ~ the coherence length, immersed 
in the superconducting (s) phase with a constant order 
parameter. This is a London-type approactP which suf- 
fices for qualitative estimates, although Ref.^argues that 
such an approach underestimates magneto-elastic effects. 

Nucleation of the normal core causes stress, since the n 
phase has a larger specific volume Vn as compared to Vs ■ 
The relative volume change C is related to the pressure 
dependence of the condensation energy or of the critical 
field 

Vn-Vs ^H, dH, 

^ Vs Att dp ' ^ ' 



The elastic energy density in isotropic solids reads.'^ 

F^\ull2 + tiul, (2) 

where Uij is the strain tensor and A, /i are Lame coef- 
ficients; summation over double indices is implied. The 
stress tensor Uij — dF/duij = XuuSij + 2fiUij, and the 
equilibrium condition daij/dxj = aijj = is given by 

Aun,j + 2^u,yj = . (3) 

For brevity, the coma in Uik.j is used to denote derivatives 
with respect to the coordinate j. 

For a single vortex directed along z, the displacement 
u = {ux,Uy,0) is radial in the plane xy, i.e., curlu — 
or u = Vx, and Uq/j = x,ai3 where x is a scalar and a, (3 
acquire only x and y values. The equilibrium condition 
([3]) reads (A -I- 2ii)x,ai3p = with the first integral 

X,i3p = V^x = C = constant . (4) 

To fix this constant, we note that x,i3i3 — Ujjp de- 
scribes compression and is related the hydrostatic pres- 
sure within the system. For the problem of the strain 
caused by a single vortex in otherwise unrestrained crys- 
tal, the pressure is zero, and we have to solve V^x = 
under the boundary condition w at large distances. 
Hence, the problem is the same as that of a linear charge 
in electrostatics: x oc Inr, r = {x,y). Hence, we obtain!^ 

rser is) ise (. 2 \ 

u, = , u^p^—^ \^6o.p - -^x^xn j ; (5) 

where is introduced for convenience and the constant 
7s is given below. 

At the core center u(0) = 0; we have 

U„ = -7nr, u'^ap = -JnSa/i (6) 

in the core interior. The constants 7's are evaluated by 
using boundary conditions at the interfaced 

^ = ^ ^ ^(^ + ^) (7) 

^" 2{x + 2^l) ' 2{x + 2^l) • ^ ' 

The displacement Ug of Eq. ^ is analogous to the elec- 
tric field of a charge with linear density 7s^^/2 situated 
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at the origin. Hence, the vortex can be considered as the 
Unear source of deformation u outside the core, whereas 
the scalar potential x satisfies 



(8) 



Vortex lattice. Consider now a 2D periodic lattice 
of vortices at positions a in an infinite sample. At first 
sight, the potential x should obey 



V2x = 2^7,e'^<5(r-a). 



(9) 



The electrostatic analogy, however, shows that this equa- 
tion cannot have a bound solution, whereas we are in- 
terested in periodic x('") to describe the infinite vortex 
lattice. We, therefore, introduce a background "charge 
density" of a sign opposite to 7sC^/2 to make the sys- 
tem "quasi-neutral" . In other words, the condition for a 
periodic x to exist is 

J dr (^27rj4^ ^ 5(r - a) C j = . (10) 

This translates to 2ttJs£.''^N + CA = where N is 
the total number of vortices and A is the area of the 
sample crossection perpendicular to the induction B; 
N/A = B/(f)Q is the density of vortices. Hence, C — 
—2Trjs(,'^B/(j)Q and we have to look for solutions of 



B 

4>o 



(11) 



an equation consistent with the equilibrium condition 
V^X = constant. 

The general solution of this equation was discussed in 
Ref.[TTl Dealing with periodic solutions, one can consider 
Xif) in a single cell under the condition dx/dn = at 
the cell boundary (n is the normal to the boundary). 
The potential within a cell centered at a = satisfies 



V^x = 2njse 



Sir) 



B 

4>o 



(12) 



The form of the unit cell depends on the vortex lat- 
tice structure which is hexagonal (triangular) in isotropic 
case of interest here. For this lattice, the boundary is a 
hexagon which - in the Wigner-Zeitz approximation - can 
be replaced with a circle. The cylindrically symmetric so- 
lution satisfying x'(-R) = with ttR^ — (j>Q/B is 



X = 7.r In 



r nB 

ro 200 



(13) 



To is an arbitrary constant irrelevant for the following. 
The crystal displacement has only one component: 



1 ttB 



(14) 



The strain tensor in cylindrical coordinated^ has two 
non-zero components: 

dur 2/^1 7r_B^ 

or \ 00 

r \r-^ 00 

The elastic energy density averaged over the cell is 
B 

Fei^— 2^r dr [\ul^{r)l2 + ^"^^(r)] , (16) 



(15) 



where the lower integration limit is the core radius on 
the order of ^. Within the London approach one cannot 
determine the radius p; we will choose it below as to 
have the elastic contribution to magnetization to be zero 
at the upper critical field Hc2- 

A straightforward evaluation gives: 



b = 



B 



c2 



and 



A = A + + p 



2^ 



(17) 



(18) 



is a quantity on the order of the elastic constants. 

Parameter C, in terms of dTc/dp. The stress de- 
pendence of the condensation energy d{H^ /8tt) /dp, to 
which the coefficient 7s is proportional, can be evaluated 
only within a detailed microscopic theory to account for 
evolution of the band structure and of the coupling re- 
sponsible for superconductivity with pressure. Such a 
calculation, if possible, would be material specific. In- 
stead, we resort to a qualitative approach to see how the 
the vortex induced strain could affect macroscopic prop- 
erties of type-II superconductors. 

First, the derivative in Eq. ([T]) can be expressed in 
terms of the measured dTc/dp: 



dp 



dp 



(19) 



Unfortunately, there is no simple enough expression for 
the condensation energy H'^/Sn = Fn — Fg for arbitrary 
temperatures, fields, and scattering regimes. The excep- 
tion is the case of a strong pair-breaking considered orig- 
inally by Abrikosov and Gor'kov,!^ who argued that due 
to extra suppression of the order parameter by, e.g., spin- 
flip scattering, the GL energy expansion holds for all tem- 
peratures down to 0. This argument has recently been 
specified for materials with zero order parameter aver- 
age (like d-wave or ±s for iron-based pnictides) .^^ This 
scheme will be used below mostly because its formal sim- 
plicity, although the qualitative results obtained have a 
broader applicability. 

The zero-field condensation energy for gapless state is 



Hi 
Stt 



A(r2 - T 



'2\2 



A. 



iV(0)r 



(20) 
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where -A''(0) is the density of states, and l/r^ = 1/t + 
I/tto whereas 1/r and l/r^ are the transport and pair- 
breaking scattering rates. One then finds: 



dTc 2ttTc 
Thus, we estimate: 

c 



t = T/T,. 



dp 



(21) 



(22) 



Also, within the gapless state, the upper critical field 
and the London penetration depth have simple T depen- 
dencies used below: 



(23) 



Magnetization. The free energy density of the mixed 
state is 



F = F^ + B''I%ti + Fl + F,, 



(24) 



where F^ is the free energy of the uniform state in zero 
field. The London energy of the vortex lattice in inter- 
mediate fields is given by 



Ft. 



In 



B 



(25) 



with 77 ^ 1. The elastic part is obtained with the help of 



Eqs. (171, (22), and (23|: 



F,^ 



B 



2'kTc dp Hc2fl 



p'^B 



c2 



(26) 



Here, A ^ 10^^ cvg/cir? is a new combination of A and fi. 

Both F]^ and fg; are evaluated here within London 
approach for Hd <^ B <^ Hc2 and, therefore, fail near 
both Hc2 and Hd- Still, we can force the magnetization 
to be zero at Hc2 by setting rj = e ~ 2.718 and p ~ 
2£_/V3. We then obtain: 



M ^ 



B-H 
An 



B OF 
4tt^ dB 

Hc2 



= Mr 



(27) 
(28) 



where is the standard London part. The average pen- 
etration length Al, that governs the field distribution in 
the mixed state, depends on the field because the aver- 
age order parameter A is suppressed by the field. This 
dependence is relevant in particular near Hc2 , since there 
the averaged order parameter ex (1 — B/Hc2) which 
translates to A cx 1/^1 - B/Hc2)-^We take this field 
dependence to characterize qualitatively Ml{T,B). 
The elastic contribution to M is 



-2A 



2TTTrH, 



c2.0 



dp 



B 1 



B 

Hc2 



(29) 



It is worth noting that this contribution is diamagnetic 
and has a minimum ai B = Hc2/2 in a reasonable agree- 
ment with data of Refs.|l]and[51 

Taking A w lO^^erg/cm^, Hc2.a = 5.9 T, Hco = 
0.35T, Al(0) w 3.5 x lO^^cm, w 10.5K, and 
dTjdp « 3 X 10^3 K cm3/erg==30 K/GPa one obtains 
M{B) shown in Fig.[l] These numbers roughly cor- 
respond to parameters for Lai.45Ndo.4oSro.i5Cu04.^ It 
should be noted that dTc/dp for this particular mate- 
rial was not measured, but the data on a similar crys- 
tal, Lai.44Ndo.4oSro.i4Cu04,'i^ show that dTc/dp\p^Q ex- 
ceeds 15K/GPa. 




FIG. 1. (Color online) The magnetization M versus applied 
field B according to Eqs. ( 27 1-( 29 1 for parameters given in 
the text. 




FIG. 2. (Color online) The magnetization M(B) measured 
in increasing and decreasing fields along the c crystal axis; 
the data are compiled from Ref.RO It is seen that M{B) is 
reversible for B > 0.4T. 

Qualitatively, the calculated M{B) is similar to 
that recorded by Ostenson et al** and shown in Fig. [2] 
The major features of the data are reproduced by the 
model remarkably well. One must bear in mind that 
the London-type isotropic model cannot pretend for 
quantitative agreement with data near Hc2 and for 
B — ^ 0. Besides, in anisotropic materials one cannot use 
the data on the dependence of hydrostatic pressure 
p as a fair representative of actual dependence of on 



4 



the stress in the plane perpendicular to the applied field. 
Moreover, the use of temperature dependencies of quan- 
tities involved characteristic of the gapless situation was 
only justified by formal simplicity of derivations. Never- 
theless, one concludes that the evidence for the vortex 
induced strain as responsible for non-monotonic behavior 
of M{B) in materials with a large dT^/dp is quite strong. 

Discussion. The elastic contribution to the vortex- 
vortex interactions was studied in a number of publica- 
tions, see Refs.[Hl IHl and references therein. It has been 
shown that this contribution is responsible for the ob- 
served flux-line lattice structures in fields tilted with re- 
spect to the c axis of NbSe2. The observed structures 
cannot be explained by London interactions alone; the 
structures seen in experiments in fact correspond to the 
maximum of the London energy. Vortex lattices are ex- 
tremely sensitive to a number of factors, among which the 
nonlocal correcti ons to London interactions were proven 
to be import antl ^^ l ^^ l Energy differences between vari- 
ous vortex lattices are exceedingly small. It is shown in 
this work that the elastic deformations caused by vortices 
may influence such a quantity as magnetization involving 
much larger energies. 

The model suggested here is profoundly qualitative. 
Materials to which the model is applied are anisotropic 
for which one needs to know a number of elastic con- 
stants. We lumped all this complexity to one number, 
A ~ 10^^ erg/cm^. Our estimates of material parame- 
ters and in particular of dTc/dp needed for evaluation 
of elastic contribution to magnetization are quite crude. 
In particular, it is hard to get reliable values of dTc/dp 
for p — >■ since commonly people are interested in high 
pressures.!^ 

We model vortices as having normal cores surrounded 
by superconductor with unperturbed order parameter so 
that the condensation energy is just H'^ /Stt which is so 



only far from _ffc2, see also Ref.'^', Hence, the model fails 
in high fields approaching Hc2- We have used simpli- 
fied T dependencies of iJc2 and corresponding to the 
strong pair-breaking situation. Still, all these uncertain- 
ties notwithstanding, the model reproduces qualitatively 
the behavior of M{B) with the second peak. 

The interpretation of the second peak in M{H) as 
an equilibrium thermodynamic property of deformable 
type-II superconductors is new; it differs from traditional 
models based on various defects-related irreversible ma- 
terial properties. The latter are always present, of course, 
and make it difficult to extract relatively weak magneto- 
elastic properties of vortex lattices. It should be noted, 
however, that well-pronounced 2nd peaks in reversible 
M[H) must result in similar peaks in irreversible mag- 
netization loops, in other words, the fish-tail loops seen 
in layered materials may also be caused by the vortex 
induced strain. 

Layered materials having strong stress dependencies of 
Tc are therefore good candidates not only for unmasking 
equilibrium magneto-elastic phenomena from the back- 
ground of strong irreversibilities, but for understanding 
the irreversible fish-tails as well. Recently, the pressure 
dependence of Tc in Ca(Fei_a:Coa:)2As2 was found to 
reach dTc/dp ~ —60 K/GPa, which is by a factor 
of 100 more than in "conventional" superconductors. 
The present work suggests that magneto-elastic effects 
should be studied in materials with large dTc/dp with 
an emphasis on macroscopic magnetization, the problem 
deserving more experimental and theoretical attention. 
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